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Abstract 

The next-to-leading order chiral pion-nucleon Lagrangian contains seven finite low- 
energy constants. Two can be fixed from the nucleon anomalous magnetic moments 
and another one from the quark mass contribution to the neutron-proton mass 
splitting. We find a set of nine observables, which to one loop order do only depend 
on the remaining four dimension two couplings. These are then determined from a 
best fit. We also show that their values can be understood in terms of resonance 
exchange related to A excitation as well as vector and scalar meson exchange. In 
particular, we discuss the role of the fictitious scalar-isoscalar meson. We also 
investigate the chiral expansion of the two P-wave scattering volumes P{~ and P^ as 
well as the isovector S-wave effective range parameter b~ . The one-loop calculation 
is in good agreement with the data. The difference P± — P^ signals chiral loop 
effects in the ttN P-waves. The calculated D- and F-wave threshold parameters 
compare well with the empirical values. 
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1 Introduction and summary 



Chiral perturbation theory is the tool to systematically investigate the consequences of 
the spontaneous and explicit chiral symmetry breaking in QCD. S-matrix elements and 
transition currents of quark operators are calculated with the help of an effective field 
theory formulated in terms of asymptotically observed fields, the Goldstone bosons and 
the low-lying baryons. A systematic perturbative expansion in terms of small external 
momenta and meson masses is possible. We call this double expansion from here on 
chiral expansion and denote the small parameters collectively by q. Beyond leading or- 
der, coupling constants not fixed by chiral symmetry appear, the so-called low-energy 
constants (LECs). For the chiral pion Lagrangian, i.e. the two-flavor case, these were 
determined more than a decade ago by Gasser and Leutwyler |I| by fitting a set of observ- 
ables calculated at next-to-leading order. In the presence of nucleons, the situation is less 
satisfactory. At next-to-leading order (q 2 ), seven finite LECs appear and 24 at order 
q 3 0, which is the first order where loops can contribute (11 of these are finite, the other 
13 are scale—dependent because they are needed to absorb the one-loop divergences). The 
dimension two pion-nucleon Lagrangian can be written as 

4 2 A = I>o,, (i) 

8=1 

with the Oi monomials in the fields of dimension two. At present, no completely systematic 
evaluation of the LECs q exists. In particular, the four LECs called 01,2,3,4 related to 
pion-nucleon scattering have been determined to one loop accuracy in the review @] 
and to order 0(q 2 ) in ||. The resulting values differ by factors of 1.5. None of these 
determinations is satisfactory since some of the input data are not very well known or 
large cancelations between individual terms appear (the best example is the isoscalar ttN 
S-wave scattering length a + ). Furthermore, if one wants to extract the dimension three 
LECs, one needs the Cj as input since they enter via 1/m suppressed vertices at that order 
(compare the form of the complete C^ N in ||). Clearly, a more stringent determination 
of these parameters is called for. A reliable determination should also be based on more 
observables than LECs in order to have some consistency checks. We close this gap in 
this paper. Without going into details, we will proceed as follows. The LECs cq and C7 
can be directly inferred from the anomalous magnetic moments of the proton and the 
neutron 0. In the absence of external (pseudo)scalar fields, the operator 5 is only non- 
vanishing for unequal light quark masses, m u 7^ m^. The corresponding LEC C5 can be 
extracted from the strong contribution to the neutron-proton mass difference. For the 
other four coupling constants, we find a set of nine observables which at one-loop order 
are given entirely in terms of tree graphs with insertions from + C^ N together with 
their 1/m-corrections and finite loop contributions. These are very special cases since in 
general at this order divergences would appear and thus dimension three LECs would be 
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needed. From a best fit to these nine observables, we are able to determine the LECs 

Cl, 2,3,4- 

Furthermore, in the meson sector it can be shown that the numerical values of the 
renormalized LECs L\(ii = M p ) can be understood to a high degree of accuracy from 
resonance saturation, i.e. they can be expressed in terms of resonance masses and coupling 
constants of the low-lying vector (V), axial-vector (A), scalar (S) and pseudoscalar (P) 
multiplets (the rf, to be precise) (in some cases, there is also some contribution from 
tensor mesons fTj). We investigate how well one can understand the numerical values of 
the Cj in terms of baryonic (A, N*, . . .) and mesonic (S, V, . . .) excitations. In particular, 
we discuss the role of the fictitious scalar-isoscalar meson and show how such correlated 
two-pion exchange reveals itself in certain LECs. Since we do not include the A as an 
active degree of freedom in the effective field theory, it contributes dominantly to some of 
the LECs as it is expected from the important role this resonance plays in pion-nuclear 
physics |§. We have already shown in a series of detailed calculations concerning a variety 
of reactions in the corresponding threshold regions that it is legitimate to encode the effects 
of the A in the pertinent LECs, see the review [HI. Since here we mostly consider threshold 
parameters (like scattering lengths and effective ranges), this procedure is expected to be 
sufficiently accurate. It remains to be proven by the authors who include the A as an 
active degree of freedom that their approach is equally precise (in the threshold region, 
of course). 

The pertinent results of this investigation can be summarized as follows: 

(i) We have determined the seven finite low-energy constants of the dimension two 
chiral pion-nucleon Lagrangian, C^n- We have found a set of nine observables that 
to one-loop order g 3 are given entirely in terms of tree graphs including insertions 
~ Ci, C2, C3, C4 and finite loop contributions, but with none from the 24 new LECs of 
C^. A best fit allows to pin down these LECs. The other three can be determined 
from the strong neutron-proton mass difference (05, which is only relevant in the 
case m u 7^ rrid) and from the anomalous magnetic moments of the proton and the 
neutron (c 6 , C7). The resulting values are listed in table 1 in section 4. 

(ii) We have shown that the empirical values of the LECs c\, . . . , C4 can be understood 
from resonance exchange. Assuming that C\ is saturated completely by scalar meson 
exchange, the values for C2, C3 and C4 can be understood from a combination of A, 
p and scalar meson exchange. It is remarkable that the scalar mass to coupling 
constant ratio Ms / y/gs needed to saturate the LEC c\ is in perfect agreement with 
typical ratios obtained in boson-exchange models of the NN force, where the <j- 
meson models the strong pionic correlations coupled to nucleons. There is, however, 
some sizebale uncertainty related to the A contribution as indicated by the ranges 
given in table 1. Concerning the LECs c§ and c-j related to k v and k s , we find that the 
isoscalar and isovector anomalous moments in the chiral limit can be well understood 
from neutral vector meson exchange. For the LEC C5, resonance saturation can not 
be used since there is no information on isospin-violating coupling constants. 
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(iii) Having established that resonance saturation can explain the LECs related to pion- 
nucleon scattering, we have considered the chiral expansion of the P-wave scattering 
volumes Pf and P 2 + to order q 3 . After renormalizing the appearing divergences, 
the chiral predictions agree at the few percent level with the empirical values. The 
largest uncertainty comes actually from the A(1232)-contribution. The difference 
Pi — P 2 + shows the relevance of chiral loops in the nN P-waves. 

(iv) The eight D and F-wave threshold parameters af± (7 = 2,3) are given to order q 3 
by lowest order tree and loop graphs only. The calculated values agree nicely with 
the empirical ones. 

This investigation is the first systematic attempt to pin down the low-energy constants 
of the chiral pion-nucleon Lagrangian. Clearly, more precise data are needed to sharpen 
the determination of the Cj. The present work, however, paves the way of fixing a subset 
of the dimension three LECs enumerated in ||. For that, a systematic study of ttN 
scattering to order g 3 should be performed. Such a study has recently been performed by 
Mojzis i. 



2 Effective Lagrangian at next— to— leading order 

(2) 

In this section, we briefly review the next-to-leading order pion-nucleon Lagrangian 
to fix our notation. We work in the path integral formulation of heavy baryon chiral 
perturbation theory which automatically obeys reparametrization invariance. All details 
are spelled out in @ or the review [|J. The pions are collected in the SU(2) matrix 
U(x) = u 2 (x) and the proton and the neutron in the iso-doublet N(x). With the 
four-velocity of the heavy nucleon fields and the covariant spin-operator a la Pauli- 
Lubanski, takes the form 



with 



N 



(v-D) 



2 m 



Kd-D- % -^{S-D,vu} 



2 m 

o2 



2 m 



+ci Tr(x+) + ( c 2 - -4- ) {y ■ u) 2 + c 3 u ■ u 

\ 8 m/ 

+ (c 4 + 4r) [S», S»}u,u v + c 5 (x+ - \^(X + )) 



4 m 



4 rri' 



(1+^)/+ +2(1+ k s )v 



N 



(2) 



X 



U, 



u^xu^ + ux^u , 



(3) 
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Here, Ffy R are the non-abelian field strength tensors of external left/right handed vector 
gauge fields and is defined analogously in terms of the isosinglet vector field vj*> 
necessary to generate the full electromagnetic current. _D M is the covariant derivative 
acting on the nucleons and, similarly, the one acting on the pions. Furthermore, x = 
2B A4 + . . . with J\4=dia.g(m u , m d ) the light quark mass matrix and B = \(Q\uu\0)\/F%, 
with F n = 92.4 MeV the pion decay constant. Some of the terms in eq.(0) receive 
1/m corrections from the expansion of the relativistic Dirac ttN Lagrangian. We have 
kept these explicitly one reason being that a phenomenological interpretation in terms of 
resonance exchange can not generate such terms. All parameters appearing are taken to 
be at their values in the chiral limit, i.e. 

Q=Q[1 + 0(m a q )} , (4) 

where m q denotes any one of the light quark masses or its average. In most cases, one 
has a = 1/2, exceptions being the anomalous isoscalar magnetic moment k, s and C5 with 
a = 1 (see below). In what follows, we can identify the nucleon mass and the axial-vector 

o o 

coupling constant with their physical values, m— m„ = 938.27 MeV and 9a= 9a — 1-26. 

(2) 

We will now be concerned with the numerical values of the LECs appearing in these 
are the q (i — 1, . . . , 5) as well as k s and k v . The machinery to do these calculations is 
spelled out in detail in [[|. 



3 Calculation of observables 

In this section, we calculate various observables to pin down the LECs q. The 01,2,3,4 
are all related to pion-nucleon threshold and subthreshold parameters and the much 
discussed pion-nucleon cx-term. We consider here only observables which to one loop 

order 0(q 3 ) are given by tree graphs including the Cj and finite loop corrections but have 

(3) 

no contribution from the 24 LECs of L,^. 

Consider first a subset of observables which depend on the LECs ci, C2 and C3. We 
introduce the small parameter fi = M n /m, i.e. the pion to nucleon mass ratio. Our 
notation concerning the ttN amplitudes and parameters is identical to the one used by 



Hohler |]I0" |. Calculation of the cr-term and the isospin-even scattering amplitude at and 
below threshold gives four relations (the one-loop contributions to the ttN scattering 
amplitude are collected in app. A), 



r2 *&nM$ 



a(0) = -4dM^ - ^f-f + 0(M*) , (5) 
< = 2 -^-(«-^) + S M(A + l)+0(Mi) , (6) 



'"01 



Fl V V 8vrF4 



4o = 



4 



with 



Note that the formula for aji was already derived in |4| for the so-called axial polarizability 
a a = 2oqi. Another relation can be derived from the isospin-even non-spin-flip scattering 
volume, Pi, 

P+ = —(l+fi)[4:a 33 + 2a 31 + 2a 13 + a u ) 

Consider also the real part of isospin-even nN forward amplitude close to threshold, 

ReT + (u) = 4vr^(a+ + b + f) + P+q 2 + 0(q 4 ) 
m 

= T + (M n ) + q 2 (3 + + 0(q 4 ) , (10) 

3 + = ——^-ReT + (u)\ M = P+ + + u)(b + + — ). (11) 

Here, to denotes the pion cms energy. The chiral expansions of the scattering length a + 
[\ \\ and of the range parameter (3 + take the form 

M 2 n 2 In 2 M 3 

T+(M.) = Ml + »)a+ = ( - Aci + 2c 2 - J£ + 2C3 ) + ^ + O(M^) , (12) 

^ = + + + - + ^(M. 2 ) • (13) 

The calculation leading to these results is somewhat tricky The tree terms are most easily 
evaluated by considering the relativistic pion-nucleon Lagrangian with the two couplings 
d 2 and d 2 , see ref.H. It leads to the forward scattering amplitude 

7T 7T 

with ujl the pion laboratory energy. Expanding in powers of 1/m gives the desired result. 

The LEC C4 appears in the chiral expansion of the isospin-odd spin-flip scattering 
volume P 2 ~ an d the subthreshold parameter 6q , 

P 2 = — (1 + fi)(a 33 - a 31 - a 13 + a u J 

= ^F 2 + T 2 ^ l + ") " 4^( 2 ^ + 3 ) + ° {Ml) ' ^ 

b 00 = ^+^)- 9 -^(l + g\) + 0(M 2 ). (16) 



5 



In the absence of a precise scheme to separate isospin-violating quark mass effects from 
the ones of virtual photons for dynamical processes, we use the information on the strong 
contribution to the neutron-proton mass difference to pin down C5, 

(m n - m p )( non ~ clm ) = 4c 5 B (m u - m d ) + 0(M*) = 4c 5 M 2 n mu ~ md + 0(M*) . (17) 

m u + md 

We assume here the standard scenario of spontaneous chiral symmetry breaking, i.e. 
B 3> F n . We remark that the generalized scenario with B ~ F n would lead to a vastly 
different value of C5. Other observables sensitive to this LEC are the ttN S-wave scattering 
lengths (taken not in the isospin limit) for processes involving at least one neutral pion 
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The anomalous magnetic moments appearing in the dimension two Lagrangian have 
been calculated in [fj, |J, 

°k s = k s + O(Ml) , (18) 
L = k v + 9 ^ + 0{MI). (19) 

47T 

These are related to the LECs c 6 and c 7 used there via 

C 6 =K V , c 7 = ^(k s - K v ) . (20) 

There are no one-loop corrections at order q 3 to k s since the spectral functions of the 
isoscalar electromagnetic form factors start at the three-pion cut, to — 9M^. 

4 Determination of the low— energy constants 

First, we must fix parameters. We use g^^ — 13.4 and qa as determined from the 
Goldberger-Treiman relation, = g^^F^jm = 1.32. We also have performed fits with 
the smaller g^^ = 13.05 and thus gA = 1-29. For the cx-term, we use <r(0) = 45 ± 
8 MeV fll3|| . The threshold and subthreshold parameters we take from P!0| , these are 
a+ = (-1.46 ± 0.10) M" 1 , bo = (10.36 ± 0.10) M" 2 , a+ = (1.14 ± 0.02) M" 3 , df = 
(1.12 ± 0.02) M" 3 , P x + = (3.01 ± 0.05) M" 3 , P 2 " = (1.00 ± 0.02) M~ 3 and b + = -(44 ± 
7) ■ 10~ 3 M^ 3 . For the isoscalar S-wave scattering length, we use a generous bound 
a + = (0±10) • 10 _3 M~ 1 since the Karlsruhe-Helsinki phase shifts fl4[] give a + = — (8.3± 



3.8) • 10 3 M n 1 where as the new PSI-ETHZ [15| value is small and positive, a + = (0 ... 4) 



10" 3 M~ x . Consequently, the value for f3 + follows to be (3 + = (2.36 ± 0.15) M~ 3 , adding 
the uncertainties of P x + , b + and a + in quadrature. The magnetic moments are known 
very precisely, for our purpose it suffices to take n v = 3.706 and k s = —0.120. Finally, 
we need a value for the strong neutron-proton mass difference. This has been evaluated 
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in great detail in [16| and we thus use (m n — m p )( non elm ) = (2.0 ± 0.3) MeV. We remind 



the reader that the photon cloud contribution as calculated via the Cottingham formula 
is about 0.8 MeV. The light quark mass ratio has been determined recently by Leutwyler 
17|, m d /m u = 1.8. 



With the error bars for the various observables as given above, we obtain as values of 
the Ci 

ci = -1.02 ± 0.06 GeV -1 , c 2 = 3.32 ± 0.03 GeV -1 , c 3 = -5.57 ± 0.05 GeV -1 . (21) 

for our central set of parameters. The uncertainties for 01,2,3 refer to the parabolic errors 
of the MINUIT fitting routine used. We remark that the fit prefers a negative value 
for a + and the a-term on the large side, a + = —10.1 ■ 10 -3 M -1 , cr(0) = 54.9 MeV. 
Clearly, the x 2 /dof of 3.03 shows that the input data are not all mutually consistent (to 
order g 3 ). Higher order corrections not yet calculated might remove these discrepancies. 
In particular, almost half of the total x 2 stems from P x , i.e. the error in c 2 and C3 is 
certainly larger than the one obtained from the fitting procedure. To get a more realistic 
estimate of the uncertainties for the various LECs, we have performed a fit were we have 
increased the uncertainties in all observables to ±15% leaving cr(0) and the range for a + 
as before. Considering the present status of the low-energy pion-nucleon scattering data 
basis, we consider such uncertainties as more realistic. For that fit, the x 2 /dof = 0.33 is 
much better and the resulting values are c\ = -0.93±0.09 GeV" 1 , c 2 = 3.34±0.18 GeV -1 , 
and C3 = —5.29 ± 0.25 GeV -1 . These we consider our central values as given in table 1 
(the uncertainties are rounded towards the larger side) together with the dimensionless 
couplings c' { = 2mCi, i = 1,...,5 (the pref actor 2m appears naturally in the heavy 
mass expansion). This fit leads to er(0) = 47.6 MeV and a + = —4.7 • 10 -3 M -1 . For 
comparison, the values determined in the review [[| based solely on the input from the 
(T-term, and a + from the Karlsruhe-Helsinki analysis, are c\ = — 0.87 ± 0.11 GeV -1 , 
c 2 = 3.34 ± 0.27 GeV -1 and c 3 = -5.25 ± 0.22 GeV" 1 . If we use the smaller value for 
g nN = 13.05 (i.e. g A = 1.29 from the GTR), we get Cl = (-1.01 ± 0.06) GeV" 1 , c 2 = 
(3.20 ± 0.03) GeV -1 and c 3 = (-5.45 ± 0.05) GeV -1 . The x 2 /dof = 3.74 is considerably 
worse. This is, however, not due to one observable but almost all of them contribute more 
to the total x 2 compared to the choice g n ^ = 13.4. Again, for the enlarged uncertainties 
one gets a substantially lower x 2 /dof = 0.30 for the values c\ = —0.91 ± 0.09 GeV -1 , 
c 2 = 3.25 ±0.18 GeV" 1 , and c 3 = -5.16 ±0.25 GeV -1 . Note that the tree level prediction 
&00 — 2m P 2 - is violated by 30%. With the inclusion of loop effects, however, a consistent 
value of C4 can be obtained from both observables. The same is true for the set of seven 
observables depending on 0^2,3. An omission of the loop corrections results in a ten times 
larger x 2 /dof. For the observables considered here, the loop effects are typically of the 
order of 30% to 50%, i.e. not small. It is also worth emphazising that we do not quote an 
uncertainty for Cq and C7 since the magnetic moments of the proton and the neutron have 
been determined with extreme precision. Notice that in [0 a somewhat larger value for 



C5 is obtained based on an leading order SU(3) estimate for (m n — m 



\non— elm 



V) 
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i 






Res 

cf CV 


Res 

Cj ranges 


1 


-0.93 ±0.10 


-1.74 ±0.19 


— 0.9* 


— 


2 


3.34 ±0.20 


6.27 ±0.38 


3.9 


2. ..4 


3 


-5.29 ±0.25 


-9.92 ±0.47 


-5.3 


-4.5... -5.3 


4 


3.63 ±0.10 


6.81 ±0.19 


3.7 


3.1... 3. 7 


5 


-0.09 ±0.01 


-0.17 ±0.02 






6 


5.83 




6.1 




7 


-2.98 




-3.0 





Table 1: Values of the LECs in GeV -1 and the dimensionless couplings c[ for i — 1, . . . , 5. 
The LECs Cqj are dimensionless. Also given are the central values (cv) and the ranges for 
the Q from resonance exchange as detailed in section 5. The * denotes an input quantity. 



5 Phenomenological interpretation of the low-energy 
constants 

In this section, we will be concerned with the phenomenological interpretation of the 
values for the LECs q. For that, guided by experience from the meson sector |J, we 
use resonance exchange. To be specific, consider an effective Lagrangian with resonances 
chirally coupled to the nucleons and pions. One can generate local pion-nucleon operators 
of higher dimension with given LECs by letting the resonance masses become very large 
with fixed ratios of coupling constants to masses. That procedure amounts to decoupling 
the resonance degrees of freedom from the effective field theory. However, the traces of 
these frozen particles are encoded in the numerical values of certain LECs. In the case at 
hand, we can have baryonic (N*) and mesonic (M) excitations, 

d= E cf + £ cf , (22) 

N*=A,R,... M=S,V,... 

where R denotes the Roper iV*(1440) resonance. We remark again that the c, are finite 
and scale-independent. 

We consider first scalar (S) meson exchange. The SU(2) Sim interaction can be written 

as 

C nS = S[c m Ti(x+) + c d Ti(uX)} ■ (23) 

From that, one easily calculates the s-channel scalar meson contribution to the invariant 
amplitude A(s, t, u) for elastic tctt scattering, 

4 16c M 2 

A s (s,t,u) = pi{M2 _ s) [2c m Ml+Us-^Ml)f+^^[^Ml+^-^Ml)} . (24) 
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Comparing with the SU(3) amplitude calculated in fll8l , we are able to relate the c m ^ to 
the c m ^ of || (setting M Sl = M Ss = M s and using the large- A c relations c m ^ = c m ^j v3 
to express the singlet couplings in terms of the octet ones), 

1 

with \c m \ =42 MeV and \cd\ =32 MeV 0. Assuming now that c\ is entirely due to scalar 
exchange, we get 

S _ 9S C m (26) 



Here, gs is the coupling constant of the scalar-isoscalar meson to the nucleons, Csn = 
—gs NN S. What this scalar-isoscalar meson is essentially doing is to mock up the strong 
pionic correlations coupled to nucleons. Such a phenomenon is also observed in the meson 
sector. The one loop description of the scalar pion form factor fails beyond energies of 
400 MeV, well below the typical scale of chiral symmetry breaking, A x ~ 1 GeV. Higher 



loop effects are needed to bring the chiral expansion in agreement with the data [19]. 
Effectively, one can simulate these higher loop effects by introducing a scalar meson with 
a mass of about 600 MeV. This is exactly the line of reasoning underlying the arguments 
used here (for a pedagogical discussion on this topic, see f2(J). It does, however, not 
mean that the range of applicability of the effective field theory is bounded by this mass 
in general. In certain channels with strong pionic correlations one simply has to work 
harder than in the channels where the pions interact weakly (as demonstrated in great 
detail in [|H|) and go beyond the one loop approximation which works well in most cases. 

r exchange, c\ = cf , we need 

180 MeV . (27) 



For ci to be completely saturated by scalar exchange, c\ = cf , we need 

Mo 



Here we made the assumption that such a scalar has the same couplings to pseudoscalars 
as the real ao(980) resonance. It is interesting to note that the effective a-meson in the 
Bonn one-boson-exchange potential [[H] with Ms = 550 MeV and g|/(47r) = 7.1 has 



Ms/y/gs = 179 MeV. This number is in stunning agreement with the the value demanded 
from scalar meson saturation of the LEC c\. With that, the scalar meson contribution to 
C3 is fixed including the sign, since c m Q > (see ref.Q), 

c f = _2^ = 2^ Cl = -1.40 GeV^ 1 . (28) 
3 MJ c m y ! 

The isovector p meson only contributes to C4. Taking a universal p-hadron coupling and 
using the KSFR relation, we find 

C4 = -j^- = 1.63 GeV -1 , (29) 
4m 
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using pz p = 6.1 ± 0.4 from the analysis of the nucleon electromagnetic form factors, the 
process NN — > Tin |22| |23] and the phenomenological one-boson-exchange potential for 
the NN interaction. 

We now turn to the baryon excitations. Here, the dominant one is the A(1232). Using 
the isobar model and the SU(4) coupling constant relation (the dependence on the off- 
shell parameter Z has already been discussed in [|J), the A contribution to the various 
LECs is readily evaluated, 

4 = _ cf = 2c} = Alg! - = 3.83 GeV- . (30) 
2 3 4 2[(m A - m) 2 - M*\ V ; 

These numbers we consider as our central values. Unfortunately, there is some sizeable 
uncertainty in these A contributions. Dropping e.g. the factor M 2 in the denominator 



of eq . (p0|) , the numerical value decreases to 2.97 GeV . Furthermore, making use of the 
Rarita-Schwinger formalism and varying the parameter Z, one can get sizeable changes 
in the A contributions ( e.g. = 1.89, = -3.03, cf = 1.42 in GeV" 1 for Z = -0.3). 
From this, we deduce the following ranges: = 1.9 . . . 3.8, = —3.8 ... — 3.0, = 
1.4... 2.0 (in GeV" 1 ). 

The Roper iV*(1440) resonance contributes only marginally, 

C R = 9 2 A rnR QQ5GeV -i 
8{m* 2 — m 2 ) 



c 



.a 

3 



9aR =-0.06 GeV' 1 , 
16(m* — m) 



of = -4^ = 0-12GeV" 1 , (31) 
8(m* — m) 

using R = 0.28 as obtained from the partial decay width T(N* — > Nrr) ~ 110 MeV 0. 

Putting pieces together, we have for C2, C3 and C4 from resonance exchange (remember 
that C\ was assumed to be saturated by scalar exchange) 



^Res 



+ c « = 3.83 + 0.05 = 3.88 



c R es = ct + cf + c? = -3.83 - 1.40 - 0.06 = -5.29 



3 

„Res 



c£ + c p 4 + c% = 1.92 + 1.63 + 0.12 = 3.67 , (32) 



with all numbers given in units of GeV -1 . Comparison with the empirical values listed in 
table 1 shows that these LECs can be understood from resonance saturation, assuming 
only that c\ is entirely given by scalar meson exchange. As argued before, the scalar 
meson parameters needed for that are in good agreement with the ones derived from 
fitting NN scattering data and deuteron properties within the framework of a one-boson- 
exchange model. We stress again that this a-meson is an effective degree of freedom which 
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parametrizes the strong -kit correlations (coupled to nucleons) in the scalar-isoscalar chan- 
nel. It should not be considered a novel degree of freedom which limits the applicability 
of the effective field theory to a lower energy scale. As pointed out before, there is some 
sizeable uncertainty related to the A contribution as indicated by the ranges for the cf es 
in table 1. It is, however, gratifying to observe that the empirical values are covered by 
the band based on the resonance exchange model. 

The LECs k s = -0.12 and k v = 5.83 can be estimated from neutral vector meson ex- 
change, in particular 

k s = Hijj , k>v = ftp . (33) 
Using e.g. the values from |22[, = —0.16 ± 0.01 and k p = 6.1 ± 0.4, we see that 



the isoscalar and isovector anomalous magnetic moments in the chiral limit can be well 
understood from u and p° meson exchange. It is amusing that the isovector pion cloud of 
the nucleon calculated to one loop allows to explain the observed difference between k p 
and k v . In strict vector meson dominance these would be equal. It is well known [|1(J that 
the low energy part of the nucleon isovector spectral functions can not be understood in 
terms of the p-resonance alone. 



6 Aspects of pion— nucleon scattering 

Having established that resonance saturation works rather well for the dimension two 
LECs, we proceed to calculate the chiral expansion of of the isovector S-wave effec- 
tive range parameter b~ and of the other two P-wave ttN scattering volumes up-to-and- 
including terms of order g 3 . Finally, we also work out the D and F-wave threshold 
parameters af ±) I = 2,3. Results for the subthreshold parameters which do not receive 
any contribution from are collected in app. B. 



6.1 Chiral expansion of P— wave scattering volumes 

We consider P 2 + , t ne isoscalar spin-flip scattering volume, and Pf related to the isovector 
spin non-flip amplitude ||, 

Pf = y(l + fl )^-2a 33 -a 31 + 2a 13 + a n ) = (-2.52 ±0.03) M" 3 , (34) 
P 2 + = y(l + //)(- 2a 33 + 2a 31 -a 13 + a n ) = (-2.74 ±0.03) M~ 3 . (35) 



Our aim is to see how well the empirical values given in eqs. (|34]j35| ) can be understood 
within chiral perturbation theory. For that, we have to account for Born terms, one 
loop graphs and insertions from (because of the crossing properties of these ampli- 
tudes). In contrast to the previous cases, the one-loop contributions are not finite and 
an appropriate renormalization has to be performed. 
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Consider first the Born terms. Including all terms, in particular the 7nrNN Weinberg 
vertex, the expansion to order q s gives 



Pf (Born) 
P+(Born) 



9 M— + - + — ) + 

2m 2 V M n m Am 2 j AmF 2 

n 2 

2m 2 



1 + 



Mtt_ 

2m 



M w m Am 2 J 



-2.29 Ml 



-2.22 Mf 3 (36) 
(37) 



where the numbers refer to our standard set of parameters (g n N — 13.4). We now turn 
to the chiral loop corrections at order g 3 . First, one has to perform the standard cou- 

o 

pling constant renormalization, 9 n N^> 9-kN- We use dimensional regularization and the 
corresponding renormalization scale A is varied between M p = 0.77 GeV and m* = 1.44 
GeV. In principle, this scale-dependence would be balanced by the contribution from 
the LECs. Since we use resonance saturation to estimate these, there remains a small 
scale-dependent reminder which can not be fixed (compare also ||). As a check on the 
one-loop calculation, one verifies that the divergence appearing in Pf is canceled by the 
local counter term 0\ + O2 and the one in P/ by the counter term O15 of ref. [[J . At the 
scale A = m, we have 



Pf (Loop) 



P+(Loop) 



Mtt 



48tt 2 P4 
(0.25 ± 0.05) M" 3 , 

24vr 2 PiU A 



2g\ + 5g 2 A + l) In ^ + t g \ + l - g \ + 1 



1 



7 
T 



(0.05 ± 0.02) m; 



-3 



(38) 
(39) 



The uncertainty stems from the variation in A as described above. The counter term 
contribution is estimated from A-resonance exchange employing the Rarita-Schwinger 
formalism, 



P 2 +(A) 



9nN 



M m 



m°m, 



m\(m& — 2m) ^ 2Z — 1 
2(toa — m) 2 4 



m + ttia) + Z 2 m& 



-0.35 M" 3 



9wN 



m A m, 



ra^ra 



2(m A 



m 



,2 4 



-0.33 M. 



-3 



(40) 
(41) 



for the off-shell parameter Z = —0.3. Using the non-relativistic isobar model and per- 
forming no chiral expansion one finds from the A(1232)-resonance, 



Pf (A) = P+(A) 



-0.58 m; 



(42) 



2m 2 [(m/\ — m) 2 — M 2 ] 

which is almost twice as large as before. Taking the average of both A(1232)-estimates 
and adding uncertainties in quadrature, the chiral predictions to 0(q 3 ) are 



p- = (-2.44 ± 0.13) M, 



-3 



p; 



-2.70 ± 0.12) M 



-3 



(43) 
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The major uncertainty comes here from the A(1232) contribution which seems hard to 
pin down accurately. Further contributions at 0(M n ) coming from the Roper resonance 
and the p-meson (as calculated in |24]| ) fall into the error band given in eq. ([43]) . We 
note that in both cases the Born terms are dominant and A exchange amounts to a 
18 and 17 % correction, respectively. The loop correction is very small for P 2 + an d 
roughly —10% for Pf. Interestingly the small difference between Pf and P 2 + stems 
mainly from the chiral loops. In tree level calculations |§ and also the Skyrme soliton 
model p5| , Pf — P 2 + = 4tt(1 + /i)(ai 3 — a 3 i) is actually zero as a consequence of SU(4:)- 



spin-flavor symmetry. This quantity therefore serves as an interesting signal for chiral 
loop effects in the ttN P-wave amplitudes. Furthermore the chiral expansion of these 
observables shows a good convergence (as expected for these particular P-waves scattering 
volumes). The chiral predictions are well within the empirical values for Pf and P 2 + , 
however the theoretical uncertainty is larger than the experimental one. We conclude that 
also these particular ttN threshold parameters can be understood within heavy baryon 
chiral perturbation theory. 

6.2 The isovector S— wave effective range parameter 

The real part of the isovector forward scattering amplitude T~ close to threshold takes a 
similar form than given in eqs.( |Il3| , |Tl~D for T + , with 



= dd^ML-.= Mi + /•>>- + i^ + Pf ' (44) 

The second term is proportional to the isovector S-wave scattering length, a - , which we 
already discussed in detail in f2"6fl . We just mention that the prediction given in that 



paper agrees well with the recent determinations from pionic atoms ||15|| . Therefore, we 
will discuss here the isovector S-wave effective range parameter b~, which is smaller and 
of opposite sign than the isoscalar one, b~ ~ — 0.3 6 + . We prefer to keep the kinematical 
factor 4ir(l + fi) and thus have to compare with the empirical value |ITJ. 



4tt(1 +n)b~ = (0.19 ±0.09) ■ Mf 3 . (45) 

The chiral expansion of this quantity takes the following form. From the tree graphs one 
finds up to order 0{M 7r ) ) 

4.(l + ,)6-(Bo r „) = -L---^ + -^_ (2 - 5A) 

= (0.57 - 0.30 - 0.02) • M~ 3 = 0.25 ■ M~ 3 , (46) 

with the contributions of the powers M" 1 ' 0,1 given separately. One sees that the truncation 
at order M° is already in agreement with the experimental value (which has quite a large 
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error bar). As a further contribution we only mention the chiral loop correction. After 
renormalization of the pion decay constant and the irN coupling constant it reads, 



47r(l + ^) b~ (Loop) 



AT, 



487r 2 Fi 



(5^-8)ln^ + ^-l 



0.04 • Ml 



(47) 



for A = m. Varying A between M p and m*, this number changes by less than 10 %. The 
loop correction is smaller than the experimental uncertainty and this presumably holds 
for all other order M n counter term contributions. We conclude, that also the value of 
the isovector S-wave effective range parameter b~ can be understood within heavy baryon 
CHPT. 



6.3 D— and F— wave threshold parameters 

Finally, we discuss the eight D- and F-wave threshold parameters af ± , I — 2, 3. To these 
only the Born (lowest order tree) graphs and two specific one loop graphs contribute 
at order q 3 , but no counter terms from . In the following expressions the loop 
contributions are the ones carrying the factor F~ 4 . We find that in most cases the chiral 
loop corrections are quite important to bring the chiral expansion close to the experimental 
values. The latter are taken from (JTTJ. We also remark that Mojzis' calculation M of these 
threshold parameters is prior to ours. The chiral expansion up to order q 3 reads: 
D-wave threshold parameters: 

-1.83- 10" 3 M- 5 

(48) 



a+ g 2 A (2 + v) 193gj 

2+ 120vrmF2M2 115200n 2 F£M n 
a+ + (exp) = (-1.8 ± 0.3) • 10~ 3 M" 5 



^< 2 + "> + 193 «1 , , = 2 .38 ■ 10^-' 



" 4807rmF2M2 11 5200tt 2 F^M^ 
. (exp) = (2.2 ± 0.3) • 10" 3 M- 5 (49) 



= AV + ri + l+MLz^l = 3. 21 . 10 -3 M - 

2+ UOnmF^M 2 14400^^714 
a 2+ (exp) = (3.2 ± 0.1) • 10" 3 M^ 5 (50) 



g 2 A (2 + a) 2 + gl(U+ 15tt) » 

+ t^ttt = -0.21 ■ 10" 3 M: 



480vrmF2M2 28800vr 3 F4M 7r 
.(exp) = (0.1 ± 0.2) • 10" 3 M- 5 (51) 
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F-wave threshold parameters: 



9 \ (^L + — — - ) = 0.29 • 1(T 3 M.- 7 



73 

" 3+ ~ UO^FfMl V^ 2 ^ 5376ttF 2 
a+ + (exp) = (0.42 ± 0.13) • lO^M" 7 (52) 

°- = mSwi {wkn - ^) = °- 06 ■ 10_3M * 7 

a+_(exp) = (0.15 ± 0.12) • lO^M" 7 (53) 



- = I / 2 + <fl(18-77T) , 2 X _ 

3+ 140ttF2MH 3360tt 2 F 2 m 2 



-3 n/r-7 

7T 



a 3+ (exp) = (-0.25 ± 0.02) • lO^M" 7 (54) 
1/4 3 + ^(27+14,U = aM . 1Q _3 M 



, . .... j -• - I I 1 I H 1 It! " ~\ Jf 7 

3 ~ ~ 840ttF 2 M3 W 840tt 2 F 2 
a -_(exp) = (0.10 ± 0.02) • 10' 3 M- 7 (55) 
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A Pion— nucleon scattering amplitude 



Here, we give explicit closed form expressions for the one-loop contribution to the ttN- 
scattering amplitude. In the center-of-mass (cms) frame the 7rA-scattering amplitude 
7T a (q) + N(p) — > n b (q') + N(p') takes the following form: 



rpba 
1 wN 



5 ba \g + (uj, t)+ia- (q' x q) h + (u, t) 



+ ie bac T c 



g-(u,t)+ia-(q'xq)h~(u,t) (A.l) 



with uj — v ■ q — v • q' the pion cms energy and t = (q — q') 2 the invariant momentum 
transfer squared. g ± (uj,t) refers to the isoscalar/isovector non-spin-flip amplitude and 
/i ± (c<j,t) to the isoscalar/isovector spin-flip amplitude. After renormalization of the pion 
decay constant F w and the pion-nucleon coupling constant g n N one finds the following 
one-loop contributions to the cms amplitudes g ± (cu,t) and /i ± (w,t) at order g 3 : 



g + (u,t) 



loop 



9 2 A 



4uP 

M 

3w 2 ' 



Ml - + {Ml - 2t) 



M n + 



2Ml 



t y/-t 

- arctan 

2M T 



(2u 2 + t - 2Ml) \{Ml - u 2 f 2 - Ml 



(A.2) 
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9 (w,t) 



loop 



^-r<3uj 2 1 -2 In— - 



— 6l}\/M% — u 2 arcsin 



CO 

MZ 



+g A (2u; 2 + t-2M 2 



5 . M OT 
m — — 

6 A 



Ml (Ml-oo 2 fl 2 . u 
-T + 3 arcsm TF 



+ 



2M 2 (l + 2g 2 A )--(l + 5g 2 A ) 



t 



AM 2 , /j4M 2 -t + y/=t 
— m 



2M n 



t 

2 ( 



1 + 5g 2 A ) In ^ + A ( 5 + i 3 <£) _ 2 M W 2 (1 + 2<&)} 



(A.4) 



, , 94 f t-AMi V-t 4qi v, „ 9N o /9 



(A.5) 



The analytic continuation above threshold u; > M n is done via the formulae 



vl — x 2 = —i\Jx 2 — 1 , arcsin x = — + % ln(x + Vx 2 — 1) . (A. 6) 

The t-dependences of the loop-amplitudes g ± (uj, t)i oop and /i ± (u;, £)i 00 p show an interest- 
ing structure, if one discards terms proportional to g\. The t-dependence of h + {oj, t)\ oop is 
then given by (2t — M%) / (3M%F 2 ) cr(t) loop , with a(t) the nucleon scalar form factor. Fur- 
thermore, the t-dependence of g~(u>, t)\ oop becomes equal to u/(2F 2 ) G%(t)i oop , with G^(t) 
the nucleon isovector electric form factor (normalized to unity). Finally, h~(u,t)i oop has 
the same t-dependence as — l/(4mF 2 ) G\j(t)\ oop , with G\j{t) the nucleon isovector mag- 
netic form factor. The one-loop calculation of these nucleon form factors can be found in 
[2]- 

B Results for some subthreshold parameters 

Here, we collect the results for those coefficients of the subthreshold expansion (around 
v — t — 0) which to order q 3 are pure loop effects. The experimental values are taken 
from [pi 
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Quantity 


One loop result 


Experimental value 


dtx 


g\/{MnF*M^) = 0.08M~ 5 


(0.17 ± 0.01) M" 5 


rl + 
"20 


(12 + 5^)/(1927rF*M ff ) = 0.235 M" 5 


(0.200 ± 0.005) M~ 5 


tt 02 


193^/(15360^]^) = 0.036 M~ 5 


(0.036 ± 0.003) M~ 5 


&0 + l 





(0.18 ± 0.01) M" 5 


bfo 


#>i/(607r 2 F 4 M 2 ) = 0.18 M- 5 


(-1.00 ± 0.02) M~ 5 


bi x 





(0.08 ± 0.01) M~ 7 


°20 


^m/^lOvr^M 4 ) = 0.05 M- 7 


(-0.31 ± 0.02) Af- 7 


°02 





-0.01 M" 7 


d n 


g\/{2AQu 2 F^Ml) = 0.007 M~ 6 


(-0.042 ± 0.003) M~ 6 


°^20 


(7 + £ 4 )/(168tt 2 F 4 M 2 ) = 0.032 M~ 6 


(—0.039 ± 0.002) M~ 6 




(1 + 7#i)/(19207r 2 F 7T 4 M 2 ) = 0.004 M" 6 


(0.010 ± 0.001) M~ 6 


box 


^m/(967ri^M w ) = 0.20 M" 4 


(0.24 ± 0.01) M" 4 


b w 


^ 4 m/(32vrF 4 M 7r ) = 1.06M" 4 


(1.08 ± 0.05) M" 4 


bn 





(-0.055 ± 0.005) M~ 6 


bw 


g\m/ (I92ir F*M*) = 0.18 M" 6 


(0.29 ± 0.02) M" 6 


bw 


g\mj (1920vr F 4 M^) = 0.010 M~ 6 


(0.025 ± 0.002) M~ 6 



Obviously, only in some cases the one-loop result is in good agreement with the em- 
pirical values as deduced from the Karlsruhe-Helsinki (KH) phase shift analysis. Note, 
however, that recent low energy 7riV-scattering data from PSI pi] show some disagreement 
with the KH80 solution of ttN dispersion analysis. It therefore seems necessary to redo 
the 7rN- dispersion analysis with the inclusion of these new data. A new determination of 
the subthreshold coefficients is now also called for. 

References 

[1] J. Gasser and H. Leutwyler, Ann. Phys. (NY) 158 (1984) 142. 

[2] V. Bernard, N. Kaiser, J. Kambor and Ulf-G. Meifiner, Nucl. Phys. B388 (1992) 
315. 

[3] G. Ecker and M. Mojzis, Phys. Lett. B365 (1996) 312. 

[4] V. Bernard, N. Kaiser and Ulf-G. Meifiner, Int. J. Mod. Phys. E4 (1995) 193. 
[5] V. Bernard, N. Kaiser and Ulf-G. Meifiner, Nucl. Phys. B457 (1995) 147. 
[6] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B321 (1989) 311. 



17 



[7] D. Toublan, Phys. Rev. D53 (1996) 6602. 

[8] T.E.O. Ericson and W. Weise, "Pions in Nuclei", Clarendon Press, Oxford, 1988. 

[9] M. Mojzis, "Elastic 7rN scattering to order 0(p 3 ) in heavy baryon chiral perturbation 
theory" , [hep-ph / 97mmnnn] . 

[10] G. Holder, in Landolt-Bornstein, Vol. 9b2, ed H. Schopper (Springer, Berlin, 1983). 

[11] V. Bernard, N. Kaiser and Ulf-G. Meifiner, Phys. Lett. B309 (1993) 421. 

[12] S. Weinberg, Trans. N.Y. Acad, of Set. 38 (1977) 185; 

S. Weinberg, in "Chiral Dynamics: Theory and Experiment", A.M. Bernstein and 
B.R. Holstein (eds.), Springer Verlag, Berlin, 1995. 

[13] J. Gasser, H. Leutwyler and M.E. Sainio, Phys. Lett. B253 (1991) 252, 260. 

[14] R. Koch, Nucl. Phys. A448 (1986) 707. 

[15] D. Chatellard et al., Phys. Rev. Lett. 74 (1995) 4157; 
D. Sigg et al, Phys. Rev. Lett. 75 (1995) 3245; 
D. Sigg et al, Nucl. Phys. A609 (1996) 269. 

[16] J. Gasser and H. Leutwyler, Phys. Rep. C87 (1982) 77. 

[17] H. Leutwyler, Nucl. Phys. B337 (1990) 108. 

[18] V. Bernard, N. Kaiser and Ulf-G. Meifiner, Nucl. Phys. B364 (1991) 283. 

[19] J. Gasser and Ulf-G. Meifiner, Nucl. Phys. B357 (1991) 90 

[20] Ulf-G. Meifiner, Comm. Nucl. Part. Phys. 20 (1991) 119. 

[21] R. Machleidt, K. Holinde and Ch. Elster, Phys. Rep. C149 (1987) 1. 

[22] P. Mergell, Ulf-G. Meifiner and D. Drechsel, Nucl. Phys. A596 (1996) 367. 

[23] G. Hdhler and E. Pietarinen, Nucl. Phys. B95 (1975) 210. 

[24] R. D. Peccei, Phys. Rev. 176 (1968) 1812. 

[25] G. Holzwarth and B. Schwesinger, Rep. Prog. Phys. 49 (1986) 825. 

[26] V. Bernard, N. Kaiser and Ulf-G. Meifiner, Phys. Rev. C52 (1995) 2185. 

[27] Ch. Joram et al., Phys. Rev. C51 (1995) 2144, 2159. 



18 



